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Abstract 


Recent results of the DAMA/Nal experiment for WIMP direct detection 
point to a possible annual modulation effect in the detection rate. We show 
that these results, when interpreted in terms of a relic neutralino, are compat¬ 
ible with supergravity models. Together with the universal SUGRA scheme, 
we also consider SUGRA models where the unification condition in the Higgs 
mass parameters at GUT scale is relaxed. 
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I. INTRODUCTION 


The new DAMA/Nal data (running period ^ 2) provide a further indication of a 
possible annual modulation effect in the rate for WIMP direct detection, already singled 
out by the same Collaboration using former data (running period # 1) @]. In Ref. @ 
we interpret the total sample of new and former data in terms of a relic neutralino in the 
framework of the Minimal Supersymmetric extension of the Standard Model (MSSM) [Q, 
by extending the analysis that we performed previously [Q about the DAMA/Nal results 
of the running period 1. The MSSM scheme represents a very versatile approach for 
discussing supersymmetric phenomenology at the Electro-Weak (EW) scale, and does not 
bear on too strong theoretical assumptions at higher energies. In [Q, using MSSM, we prove 
that the annual modulation data are quite compatible with a relic neutralino which may 
make up the major part of dark matter in the Universe and that some of the most relevant 
supersymmetric properties are explorable at accelerators in the near future. 

In the present paper we show that the supersymmetric features, implied by the 
DAMA/Nal modulation data, are also compatible with more ambitious supersymmetry 
schemes, where the previous phenomenological model is implemented in a supergravity 
(SUGRA) framework, especially if the unification conditions, which are frequently imposed 
at the Grand Unification (GUT) scale, are appropriately relaxed [Q. 


II. SUGRA MODELS 

We remind here that the essential elements of SUGRA models are: a Yang-Mills 
Lagrangian, the superpotential, which contains all the Yukawa interactions between the 
standard and supersymmetric fields, and the soft-breaking Lagrangian, which models the 
breaking of supersymmetry. Here we only recall the soft supersymmetry breaking terms 

“ ^soft= \ 

i 

+ {+ AiXtQ.HiDt + A:XiQ<.H2Ui + h.c 

i 

where the <pi are the scalar helds, the Aj are the gaugino helds. Hi and H 2 are the two Higgs 
fields, Q and L are the doublet squark and slepton fields, respectively, and f/, D and R 
denote the S'U(2)-singlet helds for the up-squarks, down-squarks and sleptons. In Eq.(2), 
rrii and Mj are the mass parameters of the scalar and gaugino helds, respectively, and A 
and B denote trilinear and bilinear supersymmetry breaking parameters, respectively. The 
Yukawa interactions are described by the parameters h, which are related to the masses 
of the standard fermions by the usual expressions, e.g. = h^V 2 and rrib = h^Vi, where 
Vi =< Hi > are the v.e.v.’s of the two Higgs helds. 

It is worth recalling that one attractive feature of the model is the connection between 
soft supersymmetry breaking and Electro-Weak Symmetry Breaking (EWSB), which would 
then be induced radiatively. 


-BfiHiH2 + h.c.} 

( 1 ) 
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It is customary to implement the supergravity framework with some restrictive assump¬ 
tions about unihcation at grand unihcation scale Mcut'- 

i) Unification of the gaugino masses: Mi^Mcut) = ^ 1 / 2 , 

ii) Universality of the scalar masses with a common mass denoted by tuq: 
mi{MGUT) = mo, 

iii) Universality of the trilinear scalar couplings: A\Mgut) = A!^[Mgut) = A'^[Mgut) 

= ^omo- 


As extensively discussed in Ref. [^, these conditions have strong consequences for low- 
energy supersymmetry phenomenology, and in particular for the properties of the neutralino, 
which is defined as the lowest-mass linear superposition of the two neutral gauginos (7 and 
Z) and the two neutral higgsinos {Hi and H 2 )'- 

X = oi 7 + CI2Z -|- (I3H1 CI4H2 ■ ( 2 ) 


Different neutralino compositions are classihed in terms of the parameter P = of + a^- 
gaugino-like when P > 0.9, mixed when 0.1 < P < 0.9 and higgsino-like when P < 0.1. 

The unihcation conditions represent a theoretically attractive possibility, which makes 
strictly universal SUGRA models very predictive. However, the above assumptions, par¬ 
ticularly ii) and iii), are not fully justihed, since universality may occur at a scale higher 
than Mgut, he. the Planck scale or string scale, in which case renormalization above Mgut 
weakens universality in the m*. Deviations from some of the unihcation conditions have 
been considered by a number of authors Ml. Implications of these deviations for relic 
neutralino phenomenology have been discussed in detail in Refs. [la¬ 


in the present paper, we discuss the DAMA/Nal data both in a SUGRA model with 
strict unihcation conditions and in a SUGRA framework, where we introduce a departure 
from universality in the scalar masses at Mgut which splits the Higgs mass parameters Mhi 
and Mh 2 in the following way: 


^hX^gut) — mo(l + di) . (3) 

The parameters di which quantify the departure from universality for the Mjj^ will be var¬ 
ied in the range (— 1 ,- 1 - 1 ), but are taken to be independent of the other supersymmetric 
parameters. 

Our supersymmetric parameter space is constrained by the following conditions: (a) 
all experimental bounds on Higgs, neutralino, chargino and sfermion masses are satished 
(for current LEP bounds see, for instance. Ref. ||T^Jl3| ), (b) the neutralino is the Lightest 
Supersymmetric Particle (LSP), (c) the constraints on the 6 —*• s -|- 7 process are satished, 
(d) the constraints on the mass of the bottom quark rub are also satished, however with a 
b — T Yukawa unihcation relaxed by about 20%, (e) EWSB is realized radiatively, (f) the 
neutralino relic abundance does not exceed the cosmological upper bound, which is 

conservatively set here as < 0.7. Because of the requirements of radiative EWSB 

and of the universality conditions on the gaugino masses and on the trilinear couplings, the 
independent supersymmetric parameters are reduced to the following set (apart from the 
5i’s): mi/2, mo, Ao, tan/? = V 2 /V 1 and sign(p). 
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The Renormalization Group Equations (RGB’s) are solved by using the 1-loop beta 
functions including the whole supersymmetric particle spectrum from the GUT scale down 
to Mz, neglecting the possible effects of intermediate thresholds. Two-loop and threshold 
effects on the running of the gauge and Yukawa couplings are known not to exceed 10% of 
the hnal result [113. While this is of crucial importance as far as gauge coupling unihcation 


is concerned [Q, it is a second-order effect on the evolution of the soft masses and then it 


is neglected here. 

In order to specify the supersymmetry phenomenology, boundary conditions for the gauge 
and Yukawa couplings have to be specihed. Low-scale values for the gauge couplings and 
for the top-quark and the tan-lepton Yukawa couplings are hxed using present experimental 
results. In particular, we assign for the top mass the value = 175 GeV. 

A few qualifications are in order here about the constraints due to the & —> s + 7 process 
and for the bottom mass. In our analysis, the inclusive decay rate BR(i? —>■ Y^y) flUH 
is calculated with corrections up to the leading order. Next-to-leading order corrections 


21 ]^^ are included only when they can be applied in a consistent way, i.e. both to stan¬ 
dard model and to susy diagrams. This criterion limits the use of next-to-leading order 
corrections to peculiar regions of the supersymmetric parameter space, where the assump¬ 
tions, under which the next-to-leading order susy corrections have been obtained, apply 
||24|| . We require that our theoretical evaluation for BR(i? —>■ Y^y) is within the range: 1.96 
xl0“^ < BR(R —>■ Y^y) < 4.32 xl0“^. This range is obtained by combining the experi¬ 
mental data of Refs. p5| , p6| at 95% G.L. and by adding a theoretical uncertainty of 25%, 
whenever the still incomplete next-to-leading order susy corrections cannot be applied. 

The supersymmetric corrections to the bottom mass include contributions from bottom- 
squark-gluino loops and from top-squark-chargino loops |27|. In the present analysis, the 


bottom mass is computed as a function of the other parameters and required to be compatible 
with the present experimental bounds. Theoretical uncertainties in the evaluation of mb arise 
from the running of the RGB’s. Since our choice is to solve RGB’s at the 1-loop level and 
without thresholds, we estimate an uncertainty of the order of 10 % in our prediction for m^. 
To take into account such an uncertainty we have chosen to weaken the bounds on m^ given 
in P by 10%. Thus we require m^ to fall into the range 2.46 GeV < mi,{Mz) < 3.42 GeV, 
at 95% G.L. As mentioned above, b — t Yukawa unihcation is relaxed by about 20%. 

The neutralino relic abundance is calculated as illustrated in [^. As already stated, 
we apply to our supersymmetric parameter space an upper bound conservatively set at the 
value < 0.7, and we consider as cosmologically interesting the range 0.01 < < 0.7. 

H 


However, we stress that, according to the most recent data and analyses 


the most 


appealing interval for the neutralino relic abundance is the narrower one: 0.02 < VL^h? < 0 . 2 . 
The local neutralino density is factorized in terms of the total local dark matter density 
Pi as p^ = ^pi- The parameter ^ is calculated according to the usual rescaling recipe |]30| : 
^ = min(r2^h^/(r2h^)min, !)• We take here (r2h^)jnin = 0.01. 

In our analysis the mi/ 2 , mo, Aq, tan/? parameters are varied in the following ranges: 
10 GeV < mi /2 < 500 GeV, mo < 1 TeV, —3 < Aq < -|-3, 1 < tan/? < 50; the 
parameter p is taken positive. We remark that the values taken here as upper limits of 
the ranges for mi/ 2,^0 are inspired by the upper bounds which may be derived for these 
quantities in SUGRA theories, when one requires that the EWSB, radiatively induced by 
the soft supersymmetry breaking, does not occur with excessive hne tuning (see Ref. [H and 
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references quoted therein). 


III. RATE FOR WIMP DIRECT DETECTION 

The indication of a possible annual modulation effect, singled out by the DAMA/Nal 
data 1^, points, at a 2-a C.L., to a very delimited region in the plane ^ where 

is the WIMP mass and is the WIMP-nucleon elastic scalar cross section. 

When the uncertainties in pi are taken into account, the original 2-a C.L. region singled 
out by DAMA/Nal data has to be enlarged into a region R, where the quantity 
falls into the following range 0: (1 — 3) x 10“® nb < ^ 3 10“® nb, in the mass 

range 30 GeV < <110 GeV. 

In Ref. P], where the relevant formulae for the evaluation of reported, it is 

also shown that, in the MSSM, wide domains of the supersymmetric parameter space provide 
values for which are within region R. Out of the two competing contri¬ 
butions to Higgs-exchange and squark-exchange processes, usually the former 

dominates over the latter one. Then let us discuss which properties of the Higgs-exchange 
amplitude are instrumental in making this contribution sizeable enough as required by the 
modulation effect: ~ 1 ^ 10”^ As is clear from the expressions given in Ref. [Q, 

the most important parameters for establishing the size of the Higgs-exchange amplitude 
are: rrih, tan/3, and the mixing angle a of the two GP-even neutral Higgs bosons {h and H). 
The largest values for the Higgs-exchange amplitude occur for the following combination of 
their respective values: small mh (because of the propagator), large values of tan/3 and of a 
(because of the structure of the Higgs-quark couplings). In SUGRA these three parameters 
are rather strongly correlated [§], so that requiring some lower bound imposes 

severe constraints on the parameter space. These properties may be suitably discussed in 
terms of the mass niA of the GP-odd neutral Higgs boson (however, we remind that this 
parameter is not free, but depends on the parameters dehning our parameter space and on 
the RGB’s evolution). 

From Fig. 1, which displays a generic scatter plot of versus in the universal 

SUGRA model (i.e. with (5/s = 0), we see that the lower bound ■C'^scaiar^ 1 ^ 
implies the upper bound rriA < 180 GeV. If we further take into account the scatter plot of 
ruA in terms of tan/3, shown in Fig. 2, we obtain in turn a lower limit on tan/3: tan/3 > 42. 

We also note that just a small relaxation in the lower bound for ^ 1 ^ 10”^ 

would also allow intermediate values of tan/3: tan/3 ~ 5 — 10. The reasons for the typical 
feature of the scatter plot of Fig. 2 for < 300 GeV are discussed in the Appendix. 

If we now require our supersymmetric conhgurations to lie inside the modulation region 
R (this set of conhgurations is denoted as set S), we obtain the scatter plot of Fig. 3, which 
proves that the annual modulation data are in fact compatible with a universal SUGRA 
scheme. Fig. 4 shows that a number of the selected supersymmetric conhgurations fall into 
the cosmological interesting range of 

The other qualihcations for the conhgurations which he inside the region R, relevant for 
searches at accelerators, concern the ranges for the h-Higgs boson mass, the neutralino mass 
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and the lightest top-squark mass, which we hnd to be: ruh < 115 GeV, 50 GeV < < 

100 GeV and 200 GeV < mi^ < 700 GeV, respectively. 

Let ns now tnrn to SUGRA schemes with deviations from nniversal scalar masses (i.e. 
Si’s 7 ^ 0). By varying the nsnal snpersymmetric parameters as before and the Si’s in the 
range — 1 < 5, < +1 we hnd the scatter plot of Fig. 5, which shows that the reqnirement 
'Cc’'sraiar°”^ > 1 X 10“® nb implies now a more relaxed npper bonnd on tua < 330 GeV. 
From Fig. 6, we see that this npper limit on tua is compatible with all valnes of tan/?. 
In Figs. 7 and 8 we notice that, as expected, onr new sample of representative points 
covers a slightly wider domain of region R, and, more signihcantly, contains new nentralino 
conhgnrations of cosmological interest. Figs. 9-11 show how the valnes for the h-Higgs 
boson mass, the nentralino mass and the lightest top-sqnark mass are distribnted in terms 
of tan /3. The ranges of these masses are similar to those already fonnd in the nniversal case, 
bnt now tan (3 extends to the interval 10 < tan (3 < 50, instead of being limited only to very 
large valnes. 

The distribntion of the valnes for the parameters dj’s, which provide snpersymmetric 
conhgnrations in agreement with the annnal modnlation data, are shown in Fig. 12. The 
pecnliar distribntion of representative points in the left-npper side of the hgnre may be 
easily nnderstood in terms of the general properties discnssed in Ref. 0. The generic trend 
displayed in this hgnre shows that small valnes of tan/? reqnire sizeable deviations from 
nniversality in the Higgs mass parameters. 


IV. CONCLUSIONS 

We have analysed the total sample of new and former DAMA/Nal data which 

provide the indication of a possible annnal modnlation ehect in the rate for WIMP direct 
detection. We have demonstrated that these experimental data, already proved to be widely 
compatible with relic nentralinos of cosmological interest in a MSSM scheme ii, are also 
compatible with a SUGRA framework. 

We have specihcally considered a snpergravity scheme with strict nnihcation conditions 
on scalar and gangino masses and on trilinear scalar conplings, and snpergravity models with 
departnres from nniversality in the scalar mass parameters of the Higgs sector. We have 
proved that in nniversal SUGRA nentralino conhgnrations with interesting cosmological 
properties may be fonnd. Even more so, in case of non-nniversal SUGRA models. 

Other main resnlts of onr analysis are the following. In the nniversal SUGRA model 
the constraints imposed by the DAMA/Nal data imply for the h-Higgs boson mass, the 
nentralino mass and the lightest top-sqnark mass, the following ranges: rrih <115 GeV, 50 
GeV < rriy < 100 GeV and 200 GeV < mr, < 700 GeV, respectively. In nniversal SUGRA 
tan/? is constrained to be large, tan/3 > 42, whereas, with departnre from nniversality in 
the scalar masses, the range for tan/? widens to 10 < tan/3 < 50. 


V. APPENDIX 

In this Appendix we show that the pecnliar behavionr of the scatter plot in Fig. 2 for 
mA < 300 GeV is indnced, in the nniversal SUGRA model, by the combined ehect of the 
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limits on the 6 —> s + 7 decay rate and on the mass of the bottom quark m?,. Let us start by 
reminding that in generic SUGRA models is a function of the other parameters through 
the radiatively induced EWSB mechanism. A useful parametrization is 

m\ = Kiml /2 + K 2 ml + K^Alml + A'4Ao"io"ii/2 - (4) 

where the coefficients AT* are only functions of tan/3 and the Si’s. Their properties are 
extensively discussed in Ref. 0 , to which we refer for details. To simplify the discussion let 
us take Aq = 0 . 

In the universal SUGRA model the coefficients Ki and K 2 are both decreasing functions 
of tan [3, being iLi ~ 3 and K 2 ~ 1 for tan (3 ~ 5-10 and both vanishingly small at large tan /3. 
Actually, while Ki is always positive, K 2 becomes negative for tan/3 ~ 50. Furthermore 
mi /2 is limited from above by the annual modulation data which imply 110 GeV and 
then mi /2 ~ 2.5 < 270 GeV (the relation mi /2 ~ 2.5 holds, since neutralino is mainly 

gaugino-like in universal SUGRA). Instead, for mo we have a lower bound which depends 
on tan/3, as displayed in Fig. A.l. This bound arises as a combined effect of the 6 —>• s + 7 
and nib constraints. Gonsequently rriA may be small only for tan (3 < 10, where may be 
arbitrarily small, or at large tan/3, where the product K 2 niQ is kept small by K 2 - 

In SUGRA models with 7 ^ 0 the coefficient K 2 may become vanishingly small also for 
intermediate values of tan (3 so that, in non-universal models, small values of niA may occur 
over the whole tan/3 range, as shown in Fig. 6 . 
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FIGURE CAPTIONS 


Figure 1 - Scatter plot of versus for a scanning of the supersymmetric 

parameter space as defined in Sect. II, in universal SUGRA. 

Figure 2 - Scatter plot of rriA versus tan (3 for a scanning of the supersymmetric pa¬ 
rameter space as dehned in Sect. II, in universal SUGRA. 

Figure 3 - Scatter plot of set S in the plane universal SUGRA. The 

dashed contour line delimits the 2-cr G.L. region, obtained by the DAMA/Nal Gollabora- 
tion, by combining together the data of the two running periods of the annual modulation 
experiment |]^. The solid contour line is obtained from the dashed line, which refers to the 
value Pi = 0.3 GeV cm“^, by accounting for the uncertainty range of as explained in 
Sect. Ill (the region delimited by the solid line is denoted as region R in the text). The 
representative points are denoted differently depending on the values of the neutralino relic 
abundance 

Figure 4 - Scatter plot of set S in the plane in universal SUGRA. Here 

neutralinos turn out to be gaugino-like only. The two vertical solid lines delimit the 
range of cosmological interest. The two dashed lines delimit the most appealing interval for 
as suggested by the most recent observational data. 

Figure 5 - Scatter plot of ^o'scaiar°’^^ versus uia for a scanning of the supersymmetric 
parameter space as dehned in Sect. II, in non-universal SUGRA. 

Figure 6 - Scatter plot of tua versus tan (3 for a scanning of the supersymmetric pa¬ 
rameter space as dehned in Sect. II, in non-universal SUGRA. 

Figure 7 - As in Fig. 3, for non-universal SUGRA. 

Figure 8 - As in Fig. 4, for non-universal SUGRA. Here dots and crosses denote 
neutralinos of diherent composition according to the classihcation in Sect. H. 

Figure 9 - Scatter plot of set S in the plane rrih - tan/?, in non-universal SUGRA. 
The hatched region on the right is excluded by theory. The hatched region on the left is 
excluded by present LEP data at ^/s = 183 GeV. The dotted and the dashed curves denote 
the reach of LEP2 at energies y/s = 192 GeV and a/s = 200 GeV, respectively. The solid 
line represents the 95% G.L. bound reachable at LEP2, in case of non discovery of a neutral 
Higgs boson. 

Figure 10 - Scatter plot of set S in the plane - tan/3, in non-universal SUGRA. 
The hatched region on the left is excluded by present LEP data. The dashed and the solid 
vertical lines denote the reach of LEP2 and TeV33, respectively. 

Figure 11 - Scatter plot of set S in the plane mi^ - tan/?, in non-universal SUGRA. 
The hatched region is excluded by LEP data (without any restriction on other masses). 

Figure 12 - Scatter plot of the values of the (5i and 82 parameters for set S. 

Figure A.l - Scatter plot in the plane mo-tan/3, for a scanning of the supersymmetric 
parameter space as dehned in Sect. H, in universal SUGRA. 
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